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Abstract 

This paper deals with a new model for clonal network dynamics. We 
describe in detail this model and derive special equations governing im- 
mune system dynamics based on the general gradient type principles that 
can be inherent to a wide class of real living objects. A special clonal 
network is modeled by two symmetric projector matrix variables simul- 
taneously taking into account both asymmetry of the interaction to each 
other and adaptation states that can be realized owing to possible idio- 
typic clonal suppresions. We perform computer simulations of the model 
dynamics for some simple cases of relatively low dimension, paying special 
attention to the dynamics of amounts of activated receptor strings within 
clonal network. 

keywords: clonal dynamics, immune system, mathematical model, gradi- 
ent dynamical system, Lyapunov function, computer simulation, complex 
system. 

1 Introduction 

The dynamics of clonal network is a very important problem for understanding 
immune systems dynamics and starting from the first mathematical model of 
Jerne JS] a lot of mathematical models are used for modeling it (See, for ex- 
ample [241 171 rKfl 1251 FHI2T ] b It is very well known that immune system contains 
many types of B-cells which can activate each other at certain conditions. If 
the receptors of such B-cells match each other the cells become antibody cells 
and they effectively recognize antigens [241 1251 l4l 121 ) . The construction of the 
cell gene model allows one to understand how the same network can support 
different independent state of the immune system. Among these states we can 
distinguish states which are quite different in their functional activity, for exam- 
ple a virgin state when the clone state is not activated, the immune state when 
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a clone proliferates antibody cells and finally - the inhibited state when prolif- 
eration terminates. Which state is reached depends on the local topology of the 
network and conditions of gene presentation With the additional assumption 
that the parameters of the model change between early and late states one can 
understand how self-nonself interaction is accomplished by the clonal network. 
In this case several original approaches has been proposed in order to explain 
characteristic features of evolution of B-cell repertoire |23 E01 ancl evolution of 
specificity in humoral immunity based on idiotypic network 

Dynamics of this clonal type network with symmetric interactions probably 
can converge |3] to steady or oscillatory state under very general assumptions 
£Q. But in reality interactions between clones are not confined to steady states 
(attractors) but also include other almost periodic and aperiodic behaviors. 
The interaction between an excitatory and inhibitory clone is clearly asym- 
metric. Assume that a self basic gene is seen by its own clone. Recall [TT] 
that adaptation can be achieved by the suppression of this idiotypic clone by 
the corresponding anti-idiotypic clone. Once the adaptation state is achieved, 
secondary presentation of the leading gene does not give rise to any network 
response, because the gene only increases the suppression on its clone. It would 
be interesting to observe if some basic modifications to the governing network 
model equations could be made that allow really memory retention ^2^]. If it 
is not the case, then not being able to retain memory of former inhibitory gene 
encounters may not be failing of the model under regard but rather a reflection 
of the intrinsic properties of clonal network. 

The model is composed [2U ESI El CHI E3 12H| of a varying number of cell 
clones of different specificities that form a clonal network. Each clone is char- 
acterized by its specific inhibitory and excitatory receptors, which are specified 
in the model by bit strings and clonal receptor string x a G R", a S N x and 
yp £ R™, (3 G My, normalized by the conditions ||x a || = 1 and ||y^|| = 1 al- 
lowing for some components to bring either positive or negative signs. Thereby 
this makes it possible to achieve suitable adaptation state through the related 
suppression of idiotypic clone by the corresponding anti-idiotypic clone, whose 
modeling within our clonal network is realized in the framework of the pos- 
tulated self-similarity of gene action and asymmetry of clone interaction. Two 
clones can interact via soluble inhibitory genes whenever their receptor shapes, 
i.e. bit strings are complementary. Cells that become activated proliferate and 
differentiate into inhibitory secreting cells. This process takes some time during 
which another free inhibitory cells form dynamic complexes not taking often 
real role within excitatory-inhibition network dynamics. 

The existence of invariant localized dynamic patterns says that our idiotypic 
clonal network possesses certain self-structured properties. Asymmetric inter- 
actions within the network determinate network's working size and connectiv- 
ity and determine the total inhibitory cells level. Concerning the connectivity 
within the network one can conclude that the network in equilibrium selects 
for growth and retaining within the network the clones with low connectivity. 
Within this framework low connectivity is not an intrinsic property of any par- 
ticular inhibitory cell but rather is determined by the random structure of a 
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clone's receptor and the shapes of the receptors on the other clones presented 
in the system. 

The incorporation of this meta-dynamics in this excitatory-inhibitory model 
one can consider as an attempt to account for the rapid enough turn over of 
clones in the clonal network. Such processes were also studied in 1271 |S1 121 
9 within a cellular automation approach. There was found that networks of 
automata can be considered as dynamical systems being the discrete equivalent 
of differential systems. They have been recently widely used in clonal nets 
and cellular automata to model complex systems such as brain. One of the 
main advantages of these networks for a biologist involved in modeling is that 
the construction of a model requires a minimal knowledge about the numerical 
values of the parameters defining a system. The differential equation systems 
describing require biological data on cell lifetimes like thresholds for activation, 
affinity constraints of the network dynamics. Most automata models do not 
require these data, since the basic assumption is that sell populations need only 
be described by a set of some discrete values, often and 1, where means 
that a populations absent, while 1 means that it is presented at a high enough 
level. The corresponding interactions among populations are represented by 
logical functions, i.e. Boolean set function, which most often are equivalent to 
threshold automata like physical spin systems 123 03 12- Such a modeling can 
be applied to a wide class of complicated neural networks taking into account 
different states of development of the corresponding cells and also to system 
evolution differential equations discussed below. These aspects of studying our 
excitatory-inhibitory clonal type network are planned to be discussed in more 
detail in another place. In this paper we will use in part an approach devised 
in for describing evolution of the introduced clonal network. 

2 Clonal network model description 

A network under consideration models a clonal dynamics exhibiting excitatory- 
inhibitory properties. It consists of interacting cell clones generated by in- 
hibitory and excitatory genes within a fixed medium. The latter will be called 
a configuration phase space, depending strongly on the nature of interaction 
between clones. An inhibitory clone population within the network can be ef- 
fectively encoded, in general, by a real (n x n) — asymmetric projector matrix of 

n(x) 

the canonical ^U] form X := E EndM. n , where biorthogonal to each 

other vectors x a ,x a E R", a E M x C l,n, are the corresponding reciprocity 
receptor strings amplitudes responsible for an inhibitory clonal topology within 
the network and the number cardj\f x — n(x) E means exactly the amount 
of activated receptors belonging to the inhibitory clone. The inhibitory clonal 
self-similarity is realized now by means of the fundamental projector property 
X ■ X = X for all whiles of time. 

Similarly, an excitatory population can be encoded, in general, by real 
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My) 

(m x m ) -asymmetric projector matrix of the canonical form Y :— Y yp®yp € 

EndW 71 , where biorthogonal to each other excitatory vectors y0,yp G K m , 
ft G M y C l,77i, are the corresponding reciprocity receptor strings amplitudes 
responsible for an excitatory clonal receptor topology and cardM y = n(y) G Z + 
means the amount of activated excitatory receptors during the network dynam- 
ics. The excitatory clonal self-similarity is realized here also by means of the 
fundamental projector property Y • Y = Y for all whiles of time. In general, 
the integer numbers n(x),n(y) G Z + can change during the network dynamics 
because of the possible full suppression of some activated clone receptor strings. 
Concerning the interaction between the clonal populations, it is described by 

n(x) n(y) 

means of a real (m x n) -matrix Z := z payp ® %a € Hom(M. n ; M. m ) 

with parameters zp a G M, a e N x C l,n, ft £ My C l,m, responsible for 
the strengths of interaction between receptors strings of two clonal populations. 
The adaptation of some inhibitory or excitatory during interaction clones is 
modeled within our clonal network by means of a possible time dependence of 
corresponding strength parameters zp a el, a 6 M x C l,n, ft G M y C l,m, 
that the matrix Z necessarily satisfies the next important self-similar clonal 
interaction properties ZX = Z = YZ for all whiles of time. 

3 Clonal network topology and dynamics 

It is natural to endow our configuration phase space Mi XtVtZ ) C (EndW 1 xEndW' 1 ) 
x HomiW" 1 ) R m ) with a reasonable Riemannian metrics by means of the follow- 
ing scalar product on its tangent space T(M) : 

(AX, Y, Z), (X, ?, Z)^ := tr(X T X) + tr(Y T Y) + tr{Z T Z) (1) 

where X, X G T(EndW l ), Y,Y G T(EndW n ) and Z, Z G T(Hom(W l : R m )) are 
arbitrary elements of the corresponding tangent spaces. Concerning the metrics 
one can construct a gradient vector field on the projector field manifold 
M(x,y,z) generated by the Lyapunov interaction function $ : M — > K whose 
variation is 

59(X, Y, Z) := tr{Dl8Z) + tr(DjSX) + tr{D T g 5Y) (2) 

for some specified matrices Df G T(EndR n ), D g G T(EndR m ) and D h G 
T(Hom(W l ;W n )), being responsible for the asymmetry of the interaction be- 
tween clonal populations. Except the Lyapunov function variation J2J it is nec- 
essary to involve into the picture the following natural clonal phase constraints: 

tr{A T {X 2 ~X)) = 0, tr(B T (Y 2 -Y)) = 0, 

tr((ZX - Z)Q T ) = 0, tr((YZ - Z),R T ) = 0, (3) 
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holding for any A £ EndW\B £ EndW 71 and Q,R £ Hom(R n ; R m ). Con- 
straints J2J can be still augmented in many special cases by the symmetry 
conditions 

tr((XT - X)P T ) = 0, tr^y -Y)S T ) =0, (4) 

holding also for arbitrary matrices P £ EndW 1 and S £ EndW" 1 . 

Below we will consider only this strongly symmetric case of our clonal net- 
work. Concerning the constraint conditions involved above the corresponding 
gradient vector field generated by the Lyapunov function variation J2J) is given 
as 

dX/dt = [[D f ,X],X] + 

+ (Z T Z + 2I)- 1 {Z T D h - Z J ZD f )(l -X) + 
+ (1 - X)(D T h Z - D f Z 1 Z)(Z 1 Z + 2I)-\ 

dY/dt = [[D g ,Y],Y} + 

+ (1 - Y)(D h Z J - D g ZZ r ){ZZi + 2I)- 1 + 
+ (21 + ZZt)- l {ZD\ - ZZ T D g )(l - Y), 

dZ/dt = -D h X - ZD f X + 

+2(Y - l)(D h X - D g Z)(Z T Z + 2I)- 1 + 

+Z(Z T Z + 2I)- 1 (Z J D h - Z T ZD f )(l - X) (5) 

n(x) 

with the canonical ^U] representations X := ^ x a ® x a £ EndW 1 , Y := 

a£Af x 

n(y) 

Ei®!/(i£ EndW", (x a ,x a >) = <W, a, a' G J\f x C l,n, {yp,yp>) = 5pp>, 
/3eA^ 

(3,(3' £ My C l,m. Here we put also Df = Dj, D g = D T g and took the matrix 
D h £ T{End(W l ;R m )) arbitrary. 

The gradient vector field JSJ) can be more specified if to take into account the 
following Lyapunov function special case taking into account the corresponding 
self-symmetry clonal declinations: 

= tr(a h (Z ~ X(3 f Z T (3 J g Y) + tr{Z^ - Y (3 g Z (3}X)a T h ), (6) 

where (X, Y, Z) £ Mf x>ViZ \ and (3j £ EndW 1 and (3 g £ EndW n are some constant 
matrices which are close to unity matrices with respect to the corresponding 
norms in EndW 1 and EndW 1 . Then owing to the definition one finds that 

£>i X) = 2(a h - P*Ya h X), 
Df = -{pf&PlYah + alYpZp}), 
Df = -(a h Xp f Z^l+(3 g Z(3}Zal), (7) 
Similarly one can take the following Lyapunov type function 

<S>^ = tr(a h j(ZT-X(3 f ZT))+tr(a h!g (ZT-ZT(3 g Y)) + 

tr((Z - Z(3jX)alj) + tr((Z - Y^ g Z)cx\ g \ (8) 
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where (X,Y,Z) G M(x tVtZ )-, (3f G EndW 1 and (3 g G EndW n are some constant 
matrices which are close to unity matrices with respect to the corresponding 
norms in EndW 1 and EndW 1 . Then owing to the definition © , one finds easily 
that 

Df = 2a h , f (I-X/3 f ) + 2(I-P g Y)a h , g , 
Df = -(/3 f Z^a h , f + al f Zpj), 

= -{piZal ig + a h , g F0 g ). (9) 

It has to be mentioned here that the Lyapunov function © models in general 
an asymmetric case of the mutual interaction between inhibitory and excita- 
tory clone populations, realizing a really observed pattern formation structure, 
when the whole system is both under an external medication and intrinsically 
activated immune state 0- 

The special analysis still must be done concerning the possible similarities 
between receptor sets of inhibitory and excitatory genes. This means that some 
mutual relationship between projector clone operators can be realized. For 
instance, their complete weak orthogonality can be realized and the following 
additional constraint 

tr{Y^Xn T ) = (10) 

for all (X, Y, Z) G M^ x y ^ and some constant matrices £,r) G Hom(W; M. m ) 
can hold. The condition i|10|) gives rise to a little complicated gradient field 
like 10 in part reflecting the mentioned above pattern isolating property of our 
clonal populations. 

It can be also interesting to analyze the dynamical system (jHJ in the case 
when either projector matrices X G EndW 1 or Y G EndW" 1 or both ones model 
clonal populations with the fixed number of the corresponding receptor strings 
activated during the interaction between them be realized that can happen when 
the clonal network is activated artificially by means of some external medication. 
This means that integers trX — n(x) G Z + or trY = n(y) G Z + persist to 
be fixed during the system evolution. Then these integers must be conserved 
quantities for all t G M involving the additional scalar constraints. The latter 
gives rise to a little modified gradient dynamical system like (0 which we do 
not write down here, being obtained the same way as before. 

4 Spectral analysis 

Consider now two respectively biorthogonal systems of vectors x a G M", a. G 
Af x C l,n, and yp G K m ,/3 G M y C l,m, being eigenvectors of the symmetric 
projector matrices X = X T G EndW 1 and Y = Y T G EndW", respectively, 
satisfying ^H] the following conditions: 

Xx a = x a , Yyp = yp (11) 
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for all a 6 Af x C l,n and /3 G 7V y C l,m. By differentiating the equalities 
Hilt with respect to the time, one gets that 

dX/dt x a + Xdx a /dt = dx a /db, dY/dt yp + Ydyp/dt = dyp/dt (12) 

for all a G Af x C 1, n and f3 G Ay C l,m. Making use of the equations JHJ) and 
1)11(1 . one has that 

[(I-X)(D f X-XD f )+K f (I-X) + (I-X)K f ]x a = (I — X)dx a /dt, 
[(I-Y)(D g Y -YD g ) + K g (I-Y) + (I-Y)K g ] = (I — Y)dyp/dt, 

where, by definition, matrices Kf G EndW 1 , K g G EndM. m and denote the 
corresponding out commutative parts of first two equations in (jSJ). As a result 
of simple computations one gets that 

(I-X)(D f X -XD f )x a + (I-X)K f x a = (I - X)dx a /dt, 
(I -Y)(D g Y -YD g )y fi + {I -Y)K g y p = (I — Y)dyp/ dt. (13) 

From 113|) one finds easily that 

dx a /dt = (D f X - XD f )x a + (I - X)K f x a + Xz { J\ (14) 
dyp/dt = (D g Y -YD g )yp + {I-Y)K gVp + Yzf 

for some vectors zi^ G E n and z < j 3 9) G R m for all a G N x C T~n and (3 G 
My C l,m. Since, in general, tr(dX/dt) ^ and tr (dY/dt) ^ 0, we deduce 
that the integers rankX and rankY are changing in time. On the other hand, 
since trX = n(x) G Z + and trY = n(y) G Z + are integers, we see that our the 
dynamical system possesses very interesting properties related with jumping of 
the integers n(x) and n(y) G Z + at some fixed whiles of time. This phenomenon 
can be interpreted naturally as a result of activation (dis-activation) of available 
receptor strings characterizing our interacting clonal populations during the 
interaction process. Returning back to equations jQ one can observe that 
vectors Za G M" and G W l must satisfy the conditions 

< z<p,x a ,>=0, <zj 9) ,^>=0, (15) 

< DfX a , X a " > + < KfX a , Xa" >= 0, (16) 

< D gV0, V/3" > + < Kgyp, yp„ >= 

for a, a' G M x C Tji, a" M x C T~n, and for /?, j3' G N y C 1, m, 0" <£ N y C 
l,m. Relationships fTK|l can be used as some criterion for the corresponding 
ranks of the projector matrices X G EndM. m and Y G EndW"\ evolving with 
respect to the dynamical system (JjJJl, to become jumped. 
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5 Clonal network simulation 



So far introduced mathematical clonal network model has been written as a 
formal system of matrix nonlinear differential equations. In order to check the 
characteristic features of this system we have performed computer simulations 
of the model for several values of n, m £ 1+ and as well as for different initial 
conditions and forms of the corresponding constant matrices entering the sys- 
tem. The model is composed in general of n(n + 1)/2 + m(m + l)/2 + nm € Z + 
scalar differential equations, i.e. of the inhibitor share equations, excitor share 
equations and equations describing interaction between these shares. The inte- 
gration of this system requires simultaneous integrations of these scalar differ- 
ential equations, each of which is a first order ODE. Let us now illustrate the 
dynamic behavior of our system by computer simulation at some simple starting 
data. We rewrite down governing equations JSJ) in a new equivalent form which 
was used for computer simulation: 

dX/dt = [[D f ,X],X} + 

+X(XZ T YZ + 2I)- 1 X(Z 1 YD h - Z J YZXD f )(I - X) + 
+(/ - X)(D\YZ - D f XZ J YZ)X(XZ T YZ + 2iy 1 X, 

dY/dt = [[D g ,Y],Y} + 

+(I - Y)(D h XZ J - D g YZXZ 1 )Y{YZXZ 1 + 2I) -1 Y + 
Y(2J + YZXZ J )- 1 Y(ZD J h - ZXZ J YD g )(I - Y), 

dZ/dt = -D h X - YZXD S X + 

2(Y - I){D h X - D g YZ)X(XZ T YZ + 2I)~ 1 X + 
+YZX(XZ T YZ + 2)- 1 X(Z T YD h - Z T YZXD f )(I - X).(17) 

Here when writing down the system l|17() we took into account constraints 
||3J| in order to get matrix paths always lying on the corresponding projector 
parts of the manifold M. We used Matlab software to perform our simulation 
of equations (|17|l . Despite the complex structure of equations describing our 
immune clonal network, its dynamical behavior is fairly interesting and good 
modeling some real ones. 

In Fig. 1 we plotted the outcome of numerical simulations for the case 
dimX = 4, dimY = 2. The plots in Fig. 1 show that the simulations recover 
perfectly the clonal network behavior predicted by both qualitative and some 
analytic considerations. For simplicity, we first consider the matrix Dh being 
zero. Fig. 1 displays a dynamics of inhibitory network X and excitatory Y at the 
conditions and parameters indicated in the capture of the figure. The behavior 
of the diagonal elements (receptors mapping antigen concentration) of the of 
the excitor are displayed by solid line and inhibitory receptors are displayed 
by dashed line. The displayed dynamics has a simple physical meaning. The 
receptor of activator network responsible for antigen state is getting depressed 
(equal to zero) due to activation of inhibitory one (Fig. lb). At certain value 
of time the state of the inhibitory network is tending to zero also. 
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Figure 1: Behavior of the diagonal elements of matrix X - dashed line, 
Y- solid line a), traces of the matrices X and Y - b). Here D g = 
( 0.5 -0.15 0.5 0.4\ 
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(0,0,0,0)' ;yi = (l,0) J ,y 2 = (0,0) J 



It should be noted that the displayed plot strongly depends on the param- 
eters of the metrices describing the network evolution. By adjusting values of 
the elements of the matrices we can get different network dynamics. This can 
correspond to real situation in living organisms when immune system extinct 
the antigens or it doesn't and there are many of natural parameters regulating 
immune system dynamics. 

In Fig. 2 we plotted the outcome of numerical simulations for the case 
dim A" = 4 = dimK As in the previous case, for simplicity, we consider the ma- 
trix Dh being zero. The numerical solution of the dynamical system obtained 
for certain D g and Df also corresponds to the case when interaction between 
clonal variables X and Y is completely determined by matrices D g and Df 
and, namely, they change the adaptive interaction of the main matrix variables 
through the matrix variable Z. Examination of the solutions obtained reveals 
that dynamics of the diagonal elements of the matrices X and Y G EndM 4 - a) 
and b) and traces of these matrices c) is essentially nonlinear. After some quasi- 
stationary dynamics of the traces we have their sharp flips to new integer values. 
During the system evolution a new increase or decrease of X -components can 
happen which may change essentially its trace dynamics. For initial conditions 
we used two dimensional subspace of the matrices X and Y. Concerning in- 
hibitor X we chose arbitrary orthogonal vectors xi t X2 and concerning activator 
vectors we chose two orthonormal vectors. These vectors form matrices X, Y as 
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Figure 2: Behavior of the diagonal elements of matrix Y - a), X - b), traces 
of the matrices X and Y - c), and plot of all solutions of matrices X, Y and 
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Figure 3: Clonal immune system dynamics for matrix Y - a), for 
matrix X - b), Trace of the matrices X and Y - c), d correspond- 



ingly. Here D 



D h = - 



/o.o 


0.2 


0.0 


0.0 


0.0 


0.0 


\0.0 


-0.3 



/ 1 0.5 1 0.4^ 

0.5 1 0.4 1 

1 0.4 2 0.5 

\0.4 1 0.5 1 / 

/ 0.3 0.6 0.2 0.9 \ 

0.5 1 0.03 0.4 

0.2 0.5 0.1 -0.5 

\-0.5 0.03 0.04 1 J 



> D f 



{0.5 0.5 0.4^ 

1 0.3 0.1 

0.5 0.3 1 0.7 

\0.4 0.1 0.7 a 



Initial conditions: 



Z{t 



Xl 



■X'3 



Xi = (0,0,0,0) T , x 2 = (0,l,0,0) rj 



vx = (i, o, o, oy , V2 = (o, i, o, oy , Vi = Vi = (o, o, o, oy . 



ii 




trace X Strace V 



trace Y Strace Y 




300 400 500 



Figure 4: The behavior of the diagonal elements of the matrix Y - a) and X - b). 
Traces of the matrices X and Y for indexes (1) and (2) - c), d) correspondingly. 
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well at some given values of coefficients z a p do the matrix Z. The specific data 
of these procedure are presented in caption for Figure 2. From the reported plot 
we can see that nondiagonal phase trajectories (Fig. 2, d) of the matrices X , Y 
and Z have practically linear behavior and are relatively stable for large values 
of time (t ~ 1000). 

In Figure 3 a) we report the numerical simulation of all components of ma- 
trices X- a) and Y- b) and corresponding trace dynamics of these matrices X- 
c) and Y- d) for the case with new initial data and all matrices being nonequal 
to zero. In particular, we took the rank of the matrix X in initial conditions 
equal to be constant 1 and that of the matrix Y equal to 2 and initial vectors 
are standard orthogonal vectors (See Fig. 3). 

The form of the matrix Dh plays essential role in interaction of the clonal 
network. From the reported plot we can see that phase trajectories of activator 
Y are stabilized at certain value of t 6 K+ and the values of Y practically 
doesn't change in time. The simulations run until the time when the system 
either destroys owing computational errors or stabilizes. 

In Figure 4 we present the numerical simulation of the diagonal components 
of matrices X- a) and Y- b) and the corresponding trace dynamics of these ma- 
trices -c) for the network dynamics stimulated by combination of two Lyapunov 
functions J2J, JBJ and J2J), ©■ In these cases new matrices in equations are 
determined as 

D g = k g D g + D^' 2 \ Df = k f D f + D f ia) , D > h = k h D h + D [ h K2) 

but equations of the clonal dynamics are persist evidently to be the same as (|17|) . 
The first upper index (1) corresponds to Lyapunov function 10 and the second 
upper index (2) corresponds to Lyapunov function JSJ. The results of computer 
simulations for first index are presented on Figure 4 a-c) and the second one on 
the Figure 4 - d). In this case the diagonal elements X - a) diagonal elements of 
matrix Y - b) and traces of these two matrices c) demonstrate a very interesting 
nonlinear dynamics. For simplicity of representation and in order to grasp some 
characteristic features of the clonal network dynamics we used practically the 
same matrices as in the previous simulations. All parameters of matrices Df, 

D g , Dh, Dy — D^p and Dg' = D g 2 ^ are presented on captions to this figure. 
With our simulations we studied the effects of the combination influence of these 
two Lyapunov potentials. 

It was established within these simulations that coefficients before the med- 
ication matrices influence sufficiently the network dynamics. It should be noted 
that the behavior of the system with Lyapunov potential combinations of these 
two potentials (0 and JJjJ, ijHJ are much more stabilized and we have many 
opportunities by taking coefficients k g ,kf,kh small enough to obtain different 
trace dynamics of the system. Analyzing behavior of the traces describing our 
immune clonal network we can conclude that small variations of the parameters 
do not change seriously the clonal dynamics. In this case there exists some 
region of attraction when the solutions persist the same form. But at some 
values of elements of matrices these solutions change drastically. This is the 
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case happening in the real immune system dynamics when inhibitory receptors 
possess many of parameters which behave precisely in order to depress antigens 
from the clonal network. 

Within our model of an immune clonal system the matrix trace trY G Z + of 
the matrix Y, being an integer, counts the number of activated excitatory recep- 
tors strings at the moment of time t G R+ during the interaction of excitatory 
clonal sub-network with inhibitory clonal sub-network, described at the same 
moment by the matrix trace trX G Z + of the matrix X, being also an integer 
and counting, respectively, the number of activated inhibitory receptors strings 
during the interaction. As one can see, during the interaction between inhibitory 
and excitatory clones at some moments of time there are switched some new 
excitatory and inhibitory receptors strings, and further at next whiles of time 
during their interactions, some of excitatory receptors strings become depressed 
and some inhibitory receptors cither appear or become dis-activated too. This 
event can be interpreted, for instance, as follows: our inhibitory clonal network 
solved its immune task to dis-activate the excitatory clonal sub-network, and 
next at some later while of time becomes idle too, leaving itself in the initial 
awaiting state. 

Thereby, one can state that the model studied in this paper possesses many 
of properties suitable for possible responses of a real clonal network. The stimu- 
lation of the interaction between excitatory and inhibitory sub-systems demon- 
strates their expectable direct self-similar and complementary behavior. This 
model is obviously not still completely satisfactory because it needs many of 
external parameters accounting for the main important features of real clonal 
network dynamics. But we believe that the work presented in this paper is an 
alternative step towards better understanding the essence and nature of immune 
network dynamics. Since the realistic immune networks do involve much more 
than several receptor strings elements, such an analysis is a new step to a com- 
pletely novel approach to understanding the functioning of real clonal immune 
networks. 
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